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LONG-TERM  GOALS 

The  long-term  goal  of  “Radiance  in  a  Dynamic  Ocean  (RaDyO):  Radianee  and  Visibility  as  Affeeted 
by  Inherent  Optical  Properties”  is  to  perform  detailed  investigations  of  inherent  optical  properties 
(lOPs)  and  the  effects  of  lOP  variability  on  underwater  radianee  and  visibility.  As  part  of  this  effort. 
Sea  Engineering,  Ine.  (SEI)  will  also  provide  all  data  from  the  Scripps  Pier  and  Santa  Barbara  Channel 
(SBC)  RaDyO  experiments  to  the  ONR  and  to  all  other  RaDyO  investigators,  with  a  detailed  report 
deseribing  the  technical  approach  and  algorithms  used  for  data  proeessing. 

OBJECTIVES 

The  primary  objectives  of  the  Radiance  in  a  Dynamic  Ocean  (RaDyO)  program  are  to: 

1)  Examine  time-dependent  oeeanic  radiance  distribution  in  relation  to  dynamic  surface  boundary 
layer  (SBL)  processes. 

2)  Construe!  a  radianee-based  SBE  model. 

3)  Validate  the  model  with  field  observations. 

4)  Investigate  the  feasibility  of  inverting  the  model  to  yield  SBE  eonditions. 


As  part  of  the  RaDyO  projeet,  SEI  has  performed  RaDyO  lOP  data  processing  and  analysis  ineluding 
data  sharing  with  other  RaDyO  Pis.  Data  post-proeessing  ineludes  applieation  of  field-ealibrations, 
eorrections  for  temperature,  salinity,  and  seattering  effects,  and  comparisons  with  similar  lOP  data 
collected  by  RaDyO  collaborators. 

lOP  data  analysis  as  related  to  the  RaDyO  program  ineludes  detailed  investigations  of  the  differences 
between  optieal  properties  and  optical  variability  measured  off  the  R/P  FEIP  and  off  the  RW  Kilo 
Moana  (KM)  during  the  RaDyO  SBC  experiment  and  the  effeets  of  optieal  variability  on  the  modeling 
of  the  underwater  radianee  distribution  and  visibility. 

Efforts  to  quantify  relationships  between  physieal  proeesses  and  image  performanee  parameters  (point 
spread  funetion,  PSE,  and  modulation  transfer  funetion,  MTE)  and  visibility  have  just  begun.  The  aim 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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of  this  work  is  to  determine  the  feasibility  of  enhancing  or  restoring  images  modulated  by  underwater 
scattering  effects. 

APPROACH 

The  measurement  approach  for  the  Scripps  Pier  and  SBC  experiments  involved  the  use  of  a  ship-based 
optical  profiler  package  with  the  following  sensors:  SeaBird  Electronics,  Inc.  conductivity- 
temperature-depth  (CTD)  (SBE49),  WET  Tabs,  Inc.  ac-s  (86-  or  87-wavelengths),  ac-9  (with  0.2  pm 
filter),  bb-9,  ECOFL-NTU,  and  Sequoia  Scientific,  Inc.  LISST-IOOX  Type  C  (Figure  I). 


Figure  1.  Photographs  of  the  optical  profiling  package  in  operation  at  Scripps  Pier  in  January  2008. 
Individual  sensors  are  labeled  in  the  photograph  on  the  left. 

Data  collection  during  the  Scripps  Pier  experiment  was  from  14  to  24  January  2008.  Our  optical 
profiler  was  operated  on  average  five  times  per  day  between  1000  and  1600  local  time  (Pacific 
Standard  Time,  PST)  except  one  period  of  night  time  sampling,  performed  in  conjunction  with  E. 
Lenain’s  (Scripps  Institution  for  Oceanography,  SIO)  “Cadillac”  data  collection  including  light-box 
images.  The  depth  of  data  collection  ranged  between  4  and  7  m,  depending  on  tidal  conditions.  The 
participants  were:  G.  Chang,  F.  Nencioli  (University  of  California,  Santa  Barbara),  F.  Spada,  and  A. 
Whitmire  (Oregon  State  University). 

SBC  data  were  collected  from  the  R/P  FLIP  between  1 1  and  20  September  2008.  Our  optical  profiler 
was  operated  from  the  face  boom  between  1 1  and  14  September  and  then  from  the  starboard  boom  on 
the  17**'  and  20*  of  September  (Figure  2).  Profile  times  (1230  and  1600  local  time)  were  coordinated 
with  WET  Labs,  Inc.  MASCOT  package  sampling  periods  (from  the  KM).  Face  boom  casts  were 
taken  to  30  m  over  a  20  minute  sampling  period.  The  optical  profiler  was  left  in-water  between  casts 
and  shallow  (2  m)  time  series  were  collected  for  10  minutes  every  hour  between  deep  casts.  Starboard 
boom  sampling  consisted  of  only  30  m  casts  at  time  periods  when  FLIP  crew  members  were  available 
to  assist  in  data  collection.  These  deep  casts  were  also  to  30  m,  over  45  minute  sampling  periods. 


2 


In  addition  to  the  optical  profiler,  an  instrumentation  package  was  mounted  on  the  hull  of  the  FLIP,  at 
30  m.  The  sensors  on  the  package  included;  WET  Labs,  Inc.  ac-s,  ECObb-3,  and  Water  Quality 
Monitor  (WQM).  The  WQM  provides  time  series  of  temperature,  salinity,  chlorophyll  fluorescence, 
turbidity,  and  dissolved  oxygen.  The  sampling  rate  was  once  per  minute  for  the  WQM  and  once  per 
hour  for  the  other  sensors.  The  SBC  participants  were;  G.  Chang,  F.  Nencioli,  and  F.  Spada. 
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Figure  2.  (A)  Map  of  the  Santa  Barbara  Channel,  CA  with  locations  of  the  R/VKilo  Moana  (blue 
cross)  and  R/P  FLIP  (red  dot)  indicated.  (B)  Photograph  of  the  FLIP  taken  on  19  September  2008 
in  the  SBC  with  the  face  boom  and  starboard  booms  indicated.  Santa  Cruz  Island  can  be  seen  in 

background.  Photo  by  Frank  Spada. 

The  data  processing  and  analysis  approach  utilized  factory-recommended  and  published  calibration 
and  correction  procedures  for  lOP  sensors  (e.g.,  Zaneveld  et  ah,  1994;  Pegau  et  ah,  1997;  Sullivan  et 
ah,  2006;  Zhang  et  ah,  2009)  and  computations  of  optical  products  following  methods  described  by 
Boss  et  al.  (2001)  and  Twardowski  et  al.  (2001).  Underwater  radiance  fields  were  computed  from 
measured  lOPs  using  the  radiative  transfer  model.  Hydrolight.  Horizontal  visibility  was  estimated 
from  the  beam  attenuation  coefficient  using  the  simplified  visibility  equation  provided  by  Zaneveld 
and  Pegau  (2003).  PSFs  and  MTFs  were  computed  using  measured  lOPs  and  empirical  or  analytical 
solutions  described  by  Duntley  (1971)  and  Voss  (1991),  and  Hou  et  al.  (2007).  Wavelet  analysis  has 
been  performed  on  image  performance  parameters  and  physical  forcing  variables  in  order  to  determine 
the  dominant  scales  of  variability  and  coherence,  and  the  periods  at  which  they  occur  (Grinsted  et  ah, 
2004).  Data  processing  and  analysis  was  performed  by  G.  Chang  and  A.  Whitmire  (SIO  Pier  LISST). 

WORK  COMPLETED 

All  data  collected  during  the  RaDyO  SIO  Pier  and  SBC  field  experiments  have  been  processed,  quality 
assured  /  quality  controlled,  and  distributed  to  RaDyO  investigators  via  UCSB’s  website 
(http;//www. opl.ucsb.edu/ucsbradyo/siopier/siopierdata.html;  SIO  Pier  data)  and  SETs  FTP  site 
(http;//ftp. seaengineering.com/ftp/RaDyO_GCS/;  login;  SEI_Guest,  password;  SEI_Guest;  SBC  data). 

Concurrent  profiles  of  optical  properties  collected  from  the  FLIP  and  the  KM  during  the  SBC 
experiment  were  analyzed  for  optical  variability.  The  effects  of  optical  variability  on  the  underwater 
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radiance  field  and  visibility  have  been  explored.  The  paper,  “Platform  effeets  on  optieal  variability 
and  the  predietion  of  underwater  visibility,”  by  G.  Chang,  M.  S.  Twardowski,  Y.  You,  M.  Moline,  P.- 
W.  Zhai,  S.  Freeman,  M.  Slivkoff,  F.  Neneioli,  and  G.  W.  Kattawar  was  published  in  Applied  Opties  in 
May  2010. 

Optieal  properties  eolleeted  at  2-m  water  depth  from  SIO  Pier  and  in  the  SBC  were  used  to  eompute 
point  spread  funetions,  modulation  transfer  funetions,  and  horizontal  visibility.  Time  series  of  the 
image  performanee  parameters  and  visibility  were  statistieally  analyzed  with  physieal  foreing 
parameters  (wind  speed,  wind  stress,  eurrent  veloeity,  temperature,  salinity,  and  density)  and  optieal 
properties  and  produets  (lOPs,  ehlorophyll  eoneentration,  slope  of  the  speetral  partieulate  beam 
attenuation  eoeffieient  and  the  real  index  of  refraetion  of  partieles;  Boss  et  ah,  2001;  Twardowski  et 
ah,  2001)  in  order  to  determine  the  faetors  infiueneing  image  modulation  and  visibility.  The  extended 
abstraet,  “Point  spread  funetions  and  visibility:  Gaining  elarity  on  image  proeessing  in  natural  waters,” 
has  been  submitted  to  the  Oeean  Opties  Conferenee  for  a  poster  presentation. 

RESULTS 

A  diurnal  wind  pattern  was  observed  during  the  two-week  SBC  experiment.  Winds  were  generally 
ealm  (<  4  m/s)  in  the  mornings,  and  inereased  to  greater  than  6  m/s,  oftentimes  reaehing  10  m/s  by 
1600  PDT.  Starting  on  15  September,  winds  greater  than  about  5  m/s  were  sustained  over  the  eourse 
of  two  days.  The  persistent  winds  resulted  in  inereased  upper  water  eolumn  mixing,  as  evideneed  by 
the  inerease  in  30  m  temperature  and  salinity  and  the  deepening  of  the  mixed  layer  depth. 

We  show  that  the  majority  of  optieal  variability  between  the  FLIP  and  the  KM  was  due  primarily  to 
platform  effeets,  speeifieally  the  break-down  of  stratifieation  from  mixing  by  the  hull  of  R/P  FLIP 
when  winds  were  less  than  5  m/s  (Figure  3;  Chang  et  ah,  2010).  Differing  vertieal  patterns  in 
hydrographie  properties  and  dissimilar  vertieal  distributions  of  water  eonstituents  (and  therefore 
optieal  properties)  observed  between  the  KM  and  FLIP  were  the  eonsequenee  of  FLIP-indueed  mixing. 
The  FLIP-effeet  was  not  a  faetor  during  higher  wind  speeds,  as  the  upper  water  eolumn  was  well- 
mixed  at  both  platform  loeations.  Note  that  the  FLIP  was  moored  during  the  Santa  Barbara  Channel 
RaDyO  experiment.  The  FLIP  in  free-drifting  mode,  where  it  is  allowed  to  ride  with  oeeanie  motion, 
may  not  be  as  strongly  affeeted  by  wind-indueed  mixing. 

The  observed  variable  vertieal  distributions  in  optieal  properties  greatly  affeeted  modeled  underwater 
radianee  distribution  and  predietions  of  horizontal  visibility  (Figure  4).  Hydrolight-modeled 
radiometrie  properties  and  AOPs  differed  by  about  50%  between  the  KM  and  FLIP  sites  during 
stratified,  low  wind  eonditions,  as  eompared  to  less  than  10%  during  well-mixed  periods.  The 
attenuation-based  model  presented  by  Zaneveld  and  Pegau  (2003)  and  a  baekward  Monte-Carlo  model 
were  used  to  prediet  the  horizontal  visibility  of  a  blaek  target  for  high  wind  and  low  wind  eonditions. 
Results  from  both  models  were  eonsistent.  Our  results  show  that  highly  variable  lOPs  observed 
between  the  FLIP  and  KM  during  low  winds  ean  affeet  predietions  of  visibility  by  up  to  57%. 
Differenees  were  less  than  15%  between  KM  and  FLIP  predieted  visibility  range  during  well-mixed 
periods. 
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Figure  3.  (A)  Time  series  of  wind  speed  collected  from  the  KM.  Crosses  and  triangles  indicate 
periods  of  KM  optical  profiler  data  collection.  Light  gray,  vertical  dashed  lines  denote  concurrent 
KM  and  FLIP  optical  profiler  data  collection.  (B-I)  Profiles  of percent  differences  between  various 
properties  measured  from  the  KM  versus  from  the  FLIP  as  a  function  of  density,  av.  The  units  for 
ag(412),  apg(412),  Cpg(650),  and  bbp(657)  are  m'\  Crosses  indicate  periods  of  data  collection  with 
wind  speeds  >  5  m/s  and  triangles  represent  wind  speeds  <  5  m/s  (except  the  last  profile). 
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Figure  4.  Left:  Vertical  profiles  of  radiometric  and  apparent  optical  properties  modeled  using 
Hydrolight  during  (A-E)  high  wind  conditions  on  13  September  2008  at  1251  PDT  and  (F-J)  low 
wind  conditions  on  20  September  2008  at  1225  PDT.  Pluses  represent  KM  data  and y  circles  denote 
FLIP  data.  Units  for  Ed  and  Eu  are  W  m'^  nm'^for  Lu  are  W  m'^  nm'^  sr^,  and  for  rrs  is  sr^  (R  is 
dimensionless).  Right:  Comparisons  between  computed  horizontal  visibility  of  a  black  target,  y, 
using  data  collected  from  the  KM  versus  from  the  ELIP.  Crosses  and  triangles  denote  data  collected 
during  high  wind  and  low  wind  conditions,  respectively.  Dark  and  light  symbols  signify  results  for 
targets  with  a  diameter  of  0.2  and  0.3  m  (radius  of  0.1  and  0.15  m),  respectively. 

The  strong  diurnal  signal  in  winds  was  also  observed  in  time  series  of  2-m  hydrographic  and  optical 
properties  and  the  imaging  performance  parameters.  Hydrographic  properties  suggest  upper  water 
column  mixing  during  the  period  of  sustained  high  wind  speeds  (15-17  September)  and  toward  the 
end  of  the  experiment  (20  September).  The  cause  of  the  latter  mixing  event  is  unknown.  Increases  in 
scattering  properties  are  observed  during  the  two  mixing  events  (15  and  20  September;  not  shown). 
Consequently,  the  imaging  performance  parameters  and  visibility  followed  a  similar  temporal  pattern 
(180°  phase)  as  the  lOPs  and  the  wind  speed  and  wind  stress  (Figure  5).  The  slope  of  the  PSF  was 
greatly  reduced  and  visibility  decreased  during  high  wind  speeds  and  mixed  periods. 
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Figure  5.  Computed  (A)  and  (C)  PSF  and  MTF;  and  time  series  of  (B)  PS F  slope  (Voss,  1991),  (D) 
MTF  (at  3  rad^  and  50  racT^),  and  (E)  horizontal  visibility  of  a  black  target  (Zaneveld  and  Pegau, 
2003)  from  2-m  optical  time  series  collected  in  the  SBC. 


Wavelet  results  indieate  that  the  slope  of  the  PSF  and  MTF  and  horizontal  visibility  are  signilieantly 
coherent  with  all  featured  properties  except  bbp(532)  (nor  salinity,  not  shown)  at  the  diurnal  frequency 
at  the  beginning  of  the  experiment  period  (Figure  6).  This  coherency  continues  through  17  September 
for  wind  speed  and  stress  (not  shown),  density,  Cpg(532),  bp(532)  (not  shown),  hbp/hp  (not  shown),  np, 
and  Y  and  through  14  September  for  apg(532)  and  Chi.  Significant  coherence  is  also  observed  at  higher 
frequencies  (0.25  to  0.5  day)  for  all  properties  except  bbp(532)  and  salinity  between  13  and  15 
September  and/or  15  and  17  September.  Coherence  phase  arrows  indicate  that  wind  speed  and  stress 
slightly  lead  the  variability  in  imaging  performance  parameters  and  visibility.  PSF,  MTF,  and 
visibility  are  positively  related  to  Gt  and  y  and  negatively  related  to  Cpg(532),  bp(532),  and  np. 

Visibility  is,  not  surprisingly,  strongly  negatively  coherent  with  Cpg(532)  at  all  time  periods  and 
frequencies  of  data  collection  (not  shown). 

The  range  of  values  for  the  imaging  performance  parameters  and  visibility  were  similar  between  SIO 
Pier  and  SBC  (Figure  7),  which  was  surprising  given  the  dynamic  forcing  at  the  shallow-water  SIO 
Pier  site.  We  expected  much  better  image  transmission  at  the  deeper  SBC  site,  however  strong  winds 
during  the  SBC  experiment  time  period  appeared  to  have  clouded  our  hypothesis.  Although  the 
scattering  coefficients  and  imaging  performance  parameters  are  inversely  correlated,  the  highly 
scattering  minerogenic  particles  observed  at  the  SIO  Pier  site  did  not  result  in  greatly  reduced  image 
transmission  or  visibility  when  using  equations  with  the  P -approximation  (e.g.,  Hou  et  ah,  2007). 
However,  MTF  increased  significantly  during  observed  rip  currents  when  evaluating  the  MTF  with 
measured  VSFs  (Figure  7E  and  G).  Visibility  ranged  between  4  and  12  m  and  PSF  slope  varied 
between  0.8  and  1.7,  similar  to  the  values  in  the  SBC  (Figure  7D-G). 
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Figure  6.  Cross  wavelet  transforms  (XWTs)  and  wavelet  coherence  (indicated  above  each  column) 
computed  for  various  parameters  (indicated  above  each  plot)  with  PSF  slope.  Frequencies  are 
shown  in  the  y-axis  (in  days)  and  the  time  period  of  measurements  is  along  the  x-axis  (date  in 
September  2008).  The  “cone  of  influence”  (COI)  is  outlined;  regions  outside  the  COI  are  ignored 
here.  Periods  of  high  wavelet  power  with  statistical  significance  are  shown  in  deep  red  and  outlined 
in  thick,  black  lines.  In  order  to  determine  statistically  significant  coherency,  both  the  XWT  and 
wavelet  coherence  figures  must  be  interpreted  together  (i.e.  where  both  XWT  and  coherence  show 
strong,  common  frequency  peaks).  Arrows  indicate  the  phase  direction  between  the  two  properties. 

Results  from  this  analysis  suggest  that  optieal  properties,  ineluding  the  imaging  performanee 
parameters  and  visibility  were  strongly  influeneed  by  meteorologieal  proeesses  during  the  SBC 
experiment.  Inereased  wind  speeds  and  stresses  resulted  in  upper  water  column  mixing,  decreased 
water  clarity,  and  reductions  in  image  transmission  and  visibility.  Optically-derived  particle 
characteristics  such  as  relative  density  and  particle  size  distribution  (PSD)  are  also  shown  to  be  related 
to  the  variability  of  imaging  performance  parameters.  A  decrease  in  PSD  and  an  increase  in  relative 
particle  density  accompanied  reduced  image  transmission  and  visibility,  likely  due  to  the  enhanced 
effects  of  multiple  scattering  for  smaller,  denser  particles. 

The  optical  properties  measured  at  SIO  Pier  appeared  to  have  been  influenced  primarily  by  advective 
and  biological  processes.  Despite  the  dynamic  and  shallow-water  environment  at  SIO  Pier,  image 
transmission  was  found  to  be  similar  to  that  in  the  SBC.  The  presence  of  high  concentrations  of 
minerogenic  particles  observed  during  rip  currents  at  the  Pier  resulted  in  decreased  image  transmission 
when  using  measured  near-forward  VSFs  for  computations  of  the  MTF. 
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Figure  7.  Time  series  of  (A)  attenuation  at  532  nm  and  chlorophyll  concentration,  (B), 
hackscattering  at  532  nm;  and  computed  (D)  PSF  slope,  (F)  MTF  (at  3  and  50  racT^),  and  (H) 
horizontal  visibility;  (C),  PSF,  (E)  and  (G)  MTF  computed  using  ^approximation  and  measured- fi, 

respectively.  All  data  were  collected  at  2-m  at  SIO  Pier. 

IMPACT/APPLICATIONS 

Our  results  have  the  potential  to  greatly  impact  the  modeling  of  underwater  radiance  distribution  for 
the  RaDyO  project.  The  hull  of  FLIP  itself  has  the  capability  to  mix  the  upper  water  column  and  affect 
SBL  processes.  Physical  measurements  conducted  by  RaDyO  Pis  off  the  FLIP  are  likely 
compromised.  Additionally,  lOPs  were  not  collected  off  the  FLIP  during  the  Hawaii  experiment.  The 
FLIP-effect  cannot  be  evaluated  for  the  Hawaii  RaDyO  field  data  set.  Model(s)  will  be  constructed 
and/or  validated  with  physical,  radiometric,  and  apparent  optical  property  (AOP)  measurements 
collected  off  the  FLIP  and  lOP  measurements  collected  off  the  KM,  which  was  several  km  away. 

The  comprehensive  optical  data  set  collected  as  part  of  the  RaDyO  program  affords  the  opportunity  to 
make  direct  computations  of  imaging  modulation  parameters  using  field-measured  data,  including 
small-angle  VSFs.  It  also  allows  the  investigation  of  physical  and  hydrographic  effects  on  image 
modulation  in  two  different  environments. 

RELATED  PROJECTS 

“Prediction  of  Optical  Variability  in  Dynamic  Nearshore  Environments,”  funded  by  the  National 
Defense  Center  of  Excellence  for  Research  in  Ocean  Sciences  (CEROS),  is  a  related  project  (Pis: 
Chang,  Jones,  Hansen,  Twardowski,  and  Barnard).  The  objective  of  this  project  is  to  develop  a  system 
for  forecasting  marine  optical  conditions  in  the  surf  zone  for  the  purpose  of  improving  Naval 
operations.  With  our  in  situ  optical  forecast  model,  the  Navy  Elect  will  be  able  to  deploy  remote 
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drifters,  combine  drifter  data  with  meteorological  and  oceanographic  data  within  our  model,  and 
predict  optical  properties  along  a  coastline  of  interest. 

Numerical  wave  and  hydrodynamic  models  have  been  developed  and  validated  with  field 
measurements  obtained  by  moored  platforms  and  optical  floats  and  drifters  in  two  surf  zone 
environments:  Santa  Cruz,  CA  and  Waimanalo,  HI.  A  third,  validation  field  experiment  will 
commence  in  early  September  2010  in  Duck,  NC. 
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